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Abstract. 

A picture for thermodynamics of the glassy state is introduced. It assumes that one 
extra parameter, the effective temperature, is needed to describe the glassy state. This 
explains the classical paradoxes concerning the Ehrenfest relations and the Prigogine- 
Defay ratio. As a second part, the approach connects the response of macroscopic 
observables to a field change with their temporal fluctuations, and with the fluctuation- 
dissipation relation, in a generalized non-equilibrium way. 
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1. Introduction 

Thermodynamics for systems far from equilibrium has long been a field of confusion. A 
typical application is window glass, which is far from equilibrium: it is an under-cooled 
liquid that, in the glass formation process, has fallen out of its meta-stable equilibrium. 

Until recently, the general consensus reached after more than half a century 
of research was: Thermodynamics does not work for glasses, because there is no 
equilibrium.This conclusion was partly based on the failure to understand the Ehrenfest 
relations and the Prigogine-Defay ratio. It should be kept in mind that, so far, the 
approaches leaned very much on equilibrium ideas. Well known works are the 1951 
Davies- Jones paper |L], the 1958 Gibbs-DiMarzio paper [Q], while a 1981 paper by 
DiMarzio has title "Equilibrium theory of glasses" and a subtitle "An equilibrium 
theory of glasses is absolutely necessary" 0. We shall stress that such approaches 
are not always applicable, due to the inherent non-equilibrium character of the 
glassy state. Actually, the above conclusion in italics is incorrect itself, the proper 
statement being that "thermostatics does not work for glasses". Let us stress that in 
literature "thermodynamics of glasses" usually refers to actual systems under laboratory 
conditions, i.e. far from equilibrium. These terminologic issues show how regretful 
it is that "thermostatics" got lost as expression for "equilibrium thermodynamics". 

Thermodynamics is the most robust field of physics. Its failure to describe the 



Thermodynamic picture of the glassy state 



2 



glassy state is quite unsatisfactory, since up to 25 decades in time can be involved. 
Naively we expect that each decade has its own dynamics, basically independent of the 
other ones. We have found support for this point in models that can be solved exactly. 
Thermodynamics then means a description under smooth non-equilibrium conditions. 

2. Thermodynamic picture for systems having an effective temperature 

A system that, after a quench to certain low temperature, slowly relaxes to equilibrium 
is characterized by the time elapsed since the quench, sometimes called "age" or "waiting 
time" . For glassy systems this is of special relevance. For experiments on spin glasses it 
is known that non-trivial cooling or heating trajectories can be described by an effective 
age [§]. Yet we do not wish to discuss spin glasses. They have an infinity of long 
time-scales, or infinite order replica symmetry breaking. 

We shall restrict to systems with one diverging time scale, having, in the mean field 
limit, one step of replica symmetry breaking. They are systems with first-order-type 
phase transitions, having discontinuous order parameter, but usually no latent heat. 

We consider transitions for glass forming liquids as well as for random magnets. The 
results map onto each other by interchanging volume V, pressure p, compressibility k = 
—d In V/dp, and expansivity a = din V/dT, by magnetization M, field H, susceptibility 
x = (l/N)dM/dH, and "magnetizability" a = (-l/N)dM/dT, respectively. 

The picture to be investigated in this work applies to systems of which the non- 
equilibrium state involves two well-separated time-scales. It can then characterized by 
three parameters, T,p and the effective temperature T e (t). In model systems to be 
discussed below, this quantity follows from analytically solving the dynamics of the 
system; in realistic (model) glasses that can be approximately described by a two- 
timescale picture it could follow from appropriate (numerical) experiments ||. For a 
set of smoothly related cooling experiments Ti(t) at pressures Pi, one may express the 
effective temperature as a continuous function: T e> i(t) — > T e (T,p). For the given set of 
experiments this sets a surface in (T,T e ,p) space. Combining with other experiments, 
such as cooling at a different rate, or first cooling, and then heating, the surface becomes 
multi- valued. For covering the whole space one needs to do many experiments, e.g., at 
different pressures and different cooling rates. The results should agree with findings 
from heating experiments and aging experiments. Thermodynamics amounts to giving 
differential relations between observables at nearby points in this space. 

Of special interest is the thermodynamics of a thermal body at temperature T 2 in 
a heat bath at temperature Ti = T. A basic assumption is separation of time-scales, 
and consequent separation of phase space, allowing to identify entropies S\ and S2. (For 
the energy such a decomposition cannot be made.) This setup allows to maintain the 
difference in temperatures, and applies to mundane situations such as a cup of coffee, or 
an ice-cream, in a room. The change in heat of such systems obeys &Q < TidSi +T 2 dS , 2. 

A similar two-temperature approach proves to be relevant for glassy systems. The 
known exact results on the thermodynamics of systems can be summarized by the very 
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same change in heat || 0. 

dQ = TdS cp + T e dX (1) 

where S ep is the entropy of the equilibrium processes, i.e. the fast or /^-processes that 
have a timescale less than the observation time. X is the configurational entropy of the 
slow or configurational processes (a-processes), also known as information entropy or 
complexity. In the standard definition || the configurational entropy S c is the entropy 
of the glass minus the one of the vibrational modes of the crystal. For polymers this still 
includes short-distance rearrangements, which is a relatively fast mode. It was confirmed 
numerically that S c indeed does not vanish at any temperature, and it does not fit well 
to the Adam-Gibbs || relation r eq ~ exp(C /TS C ) ||. Our quantity X only involves 
long-time processes; the relatively fast ones are counted in S ep . To stress that in X only 
the slow modes of S c contribute, it would deserve a separate name; the most natural one 
being complexity. The applicability of an Adam-Gibbs- type relation r eq ~ exp(C/TX) 
remains an open issue, but it is satisfied in a toy model for the standard folklore of the 



glassy state |jTU] 



If a system has a set of processes % with very different timescales t\ and partial 
entropies Si, one can define X as the sum of the Si having Tj > r eq /10, where r eq = max r, ; 
is the equilibrium relaxation time. For a system quenched to a low temperature and 
aging there during a time t, the sum would be restricted to r, > t/10. 

2.1. First and second law 

For a glass forming liquid the first law dll = dQ + dW becomes 

dU = TdS cp + T e dl - pdV (2) 
It is appropriate to define the generalized free enthalpy 

G = U-TS cp -T e l + pV (3) 
This is not the standard free enthalpy, since T e ^ T. It satisfies 

dG= - S ep dT - XdT e + Vdp (4) 
The total entropy is 

S = S ep +X (5) 

The second law requires dQ < TdS, leading to (T e — T)dX < 0, which merely says that 
heat goes from high to low temperatures. 

Since T e = T e (T,p), and both entropies are functions of T, T e and p, the expression 

. In the glass transition region all factors, 

v 

eads to 

C P = C 1 + C 2 ^\ p (6) 

Precisely this form has been assumed half a century ago by Tool [JTI|] as starting point 
for the study of caloric behavior in the glass formation region, and has often been used 
for the explanation of experiments Uflll]. It is a direct consequence of eq. ([!]). 



([!]) yields the specific heat C p = dQ/dT 
except drT e , are basically constant. This 
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For magnetic systems the first law brings 

dU = TdS cp + T e dl - MdH (7) 

As above, one can define the generalized free energy F = U — TS cp — T e l. It satisfies 
the relation dF = —S ep dT — ldT e — MdH. In an aging system (fixed T and H) the rate 
of change F = —lT e is usually positive. In literature [§] |24| one has also considered 
the "experimental" or "dynamical" free energy F dyn = U — T(S ep +T). It evolves as 
F&yn = (T e — T)X, thus it is related to entropy production, and always negative. 



2.2. Modified Maxwell relation 



For a smooth sequence of cooling procedures of a glassy liquid, eq. (PJ) implies a 
modified Maxwell relation between macroscopic observables such as U(t,p) — > U(T,p) = 
U(T,T e (T,p),p) and V. This solely occurs since T e is a non-trivial function of T and p 
for the smooth set of experiments under consideration. 
For glass forming liquids it reads 
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T, so the right hand side vanishes. For a glassy magnet one has 
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(9) 



2.3. Ehrenfest relations and Prigogine-Defay ratio 

In the glass transition region a glass forming liquid exhibits smeared jumps in the specific 
heat C p , the expansivity a and the compressibility k. If one forgets about the smearing, 
one may consider them as true discontinuities, yielding an analogy with continuous 
phase transitions of the classical type. 

Denoting the discontinuities as AO = Oi iqui d — O g i ass , we may follow Ehrenfest and 
take the derivative of AV(T, p g (T)) = 0. This yields the "first Ehrenfest relation" 



Act 



dT 



(10) 



while for a glassy magnet Aa = AxdH g /dT. 

The conclusion drawn from half a century of research on glass forming liquids is 
that this relation is never satisfied [l]|13][14] |TjJ. This has very much hindered progress 
on a thermodynamical approach. However, from a theoretical viewpoint it is hard to 



imagine that something could go wrong when just taking a derivative. McKenna [IB 
and, independently, also we 0, have pointed out that this relation is indeed satisfied 
automatically, but it is important say what is meant by k in the glassy state. 

Let us make an analogy with spin glasses. In mean field theory they have infinite 
order replica symmetry breaking. From the early measurements of Canella and Mydosh 



17]] on AuFe it is known that the susceptibility depends logarithmically on the frequency, 



so on the time scale. The short-time value, called Zero-Field-Cooled (ZFC) susceptibility 
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is a lower bound, while the long time value, called Field-Cooled (FC) susceptibility is 
an upper bound. Let us use the term "glassy magnets" for spin glasses with one step of 
replica symmetry breaking. They are relevant for comparison with glass forming liquids. 
For them the situation is worse, as the ZFC value is discontinuous immediately below 
T g . This explains why already directly below the glass transition different measurements 
yield different values for k. These notions are displayed in figure 1. 

Previous claims about the violation of the first Ehrenfest relation can be traced 
back to the equilibrium idea that there is one, ideal k, to be inserted in fllPl) . Indeed, 
investigators always considered cooling curves V(T,pi) at a set of pressures Pi to 
determine Act and dp g /dT. However, An was always determined in another way, like 
measurement of the speed of sound, or by making pressure steps [0. In equilibrium 
such alternative determinations would yield the same outcome. In glasses this is not 
the case: the speed of sound is a short-time process, and additional pressure steps 
modify the glassy state. Therefore alternative procedures should be avoided, and only 
the cooling curves V(T,pi) should be used. They constitute a liquid surface Vji qu id(T,p) 
and a glass surface V^i ass (T, p) in (T, p, V) space. These surfaces intersect, and the first 
Ehrenfest relation is no more than a mathematical identity about the intersection line 
of these surfaces. It is therefore automatically satisfied ||. The most careful data we 
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Figure 1. Schematic plot of the field-cooled (FC) and zero-field-cooled (ZFC) 
susceptibility in realistic spin glasses and in glassy magnets, as function of temperature, 
in arbitrary units. In realistic spin glasses the infinite time or field-cooled susceptibility 
is larger than the short time or zero-field-cooled susceptibility. In magnetic analogs of 
realistic glasses the short time susceptibility even has a smeared discontinuity at the 
glass transition, yielding a value of \ that depends on the precise type of experiment 
which is performed. In glass forming liquids the same happens for the compressibility. 
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came across were collected by Rehage and Oels for atactic polystyrene [fjl| . In figure 2 
we present those data in a 3-d plot, underlining our point of view. 

The second Ehrenfest relation follows from differentiating AU(T,p g (T)) = 0. The 
obtained relation will also be satisfied automatically. However, one then eliminates 
dU/dp by means of the Maxwell relation. We have already discussed that outside 
equilibrium it is modified. The equality T e (T,p g (T)) = T implies 



dT R «9T fi 



dT P 



dp 
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Using eq. 



dT dT 
and inserting this relation, we obtain 



T g V 




dT 




(12) 



The dT/dT term constitutes the total derivative along the glass transition line, 
dT/dT + (dT/dp)dp g /dT. Its prefactor only vanishes at equilibrium (T e = T), in which 
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Figure 2. Data of the glass transition for cooling atactic polystyrene at rate 18 K/h, 
scanned from the paper of Rehage and Oels (1976): specific volume V (cm? / g) versus 
temperature T (K) at various pressures p (kbar). As confirmed by a polynomial fit, 
the data in the liquid essentially lie on a smooth surface, and so do the data in the 
glass. The first Ehrenfest relation describes no more than the intersection of these 
surfaces, and is therefore automatically satisfied. The values for the compressibility 
derived in this manner will generally differ from results obtained via other procedures. 
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case the standard Ehrenfest relation is recovered. For glassy magnets one has similarly 



AC 




dT P 



NT " dT N \~ dT 
The Prigogine-Defay ratio is defined by 

n 




AC p Ak 



(13) 



(14) 



TV(Aa) 2 

For equilibrium transitions it should be equal to unity. Assuming that at the glass 
transition a number of unspecified parameters undergo a phase transition, Davies and 
Jones derived that II > 1 H, while DiMarzio showed that in that case the correct 
value is II 



1 IS 



In glasses typical experimental values are reported in the range 
2 < II < 5. It was therefore generally expected that IT > 1 is a strict inequality arising 
from the requirement of mechanical stability. 



Since the first Ehrenfest relation is satisfied, it holds that 



n 



AO, 



1 + 



1 



1 - 



dTf, 




(15) 



TVAa(dp g /dT) VAa \ dT 

AC P and Aa can be measured by cooling at a fixed pressure, but dp g /dT depends on 
cooling experiments at two pressures, or, more precisely, on the smooth set of cooling 
experiments. Therefore dp g /dT can be small or large, implying that II depends on 
the set of experiments. As a result, it can also be below unity. Rehage-Oels found 
II = 1.09 ~ 1 at p — 1 kbar, using a short-time value for k. Reanalyzing their data 
we find from ([0|), where the proper long-time k has been inserted, a value II = 0.77. 
Notice that it is below unity. The commonly accepted inequality II > 1 is based on the 
equilibrium assumption of a unique k. Our theoretical arguments and the Rehage-Oels 
data show that this assumption is incorrect. 



2. 4- Fluctuation formula 

The basic result of statistical physics is that it relates fluctuations in macroscopic 
variables to response of their averages to changes in external field or temperature. We 
have wondered whether such relations generalize to the glassy state. 

In cooling experiments at fixed field it holds that M =M(T(t),T e (t, H), H). For a 
thermodynamic description one eliminates time, implying M —M(T, T e (T, H), H). One 
may then expect three terms: 
1 dM 



X 



x mct (t) + x loss (t) + x cont (t) 



(16) 



t:i\ 



. To find them separately, we switch 



N dH 

The first two terms add up to (l/N)(dM/dH) 
from a cooling experiment to an aging experiment at the considered T, T e and H, by 
keeping, in Gedanken, T fixed from then on. The system will continue to age, expressed 
by T e = T e (t; T, H). We may then use the equality 



dM 



dH 
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T,T e + ~df~ p 
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T,HdH 



;i7) 
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We have conjectured pO] that the left hand side may be written as the sum of fluctuation 
terms for fast and slow processes, defining the fluctuation contribution 

dM 

T,t 



x mct (t) 



OH 



(SM 2 (t)) {ast (5M 2 (t)) slow 



NT(t) NT e (t) v ' 

The first term is just the standard equilibrium expression for the fast equilibrium 
processes. Notice that the slow fluctuations enter with their own temperature, the 
effective temperature. This decomposition is confirmed by use of the fluctuation- 
dissipation relation in the form to be discussed below. From eq. (|T7l) then follows 
the "loss" term 



X loss (t) = - 



1 dM 



NdT e 



dT e 



T,t 



(19) 



It is related to an aging experiment. In some models it is small [2C][ZI], but in another 
model |10| it is of order unity. There occurs in eq. (fl6|) also a configurational term 
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(20) 



It originates from the difference in the system's structure for cooling experiments at 
nearby fields. This is the term that is responsible for the discontinuity of \ or K a ^ the 
glass transition. Its existence was anticipated by Goldstein and Jackie |L3 
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2.5. Fluctuation- dissipation relation 

Nowadays quite some attention is payed to the fluctuation-dissipation relation in the 
aging regime of glassy systems. It was put forward in works by Sompolinsky [^2] and 
Horner [ 23] |24 ] , and generalized by Cugliandolo and Kurchan |25| , see |2(J for a review. 



In the aging regime there holds a fluctuation-dissipation relation between the 
correlation function C(t, t')—< 5M{t)5M{t') > and G(t,t'), the response of <M{t)> to 
a short, small field change 5H(t') applied at an earlier time t', 
dC{t,t') 



dt' 



T e {t,t')G(t,t') 



(21) 



with T e (t,t') being an effective temperature, also denoted as T/X(t,t') [j26f . 

We have observed that in simple models without fast processes T e (t,f) = T e (t') is 
a function of one of the times only |20| |21|. One then expects that T e (t) is close to the 
"thermodynamic" effective temperature T e (t). We have shown that [^1|] 



f e (t)=T e (t)-f e (t) 



'dlnC(M') 



dt' 



t'=t 



+ 



(22) 



So the effective temperatures T e and T e are not identical. However, in the models 
analyzed so far, the difference is subleading in 1/lnt. 

Notice that the ratio d t iC(t,t')/G(t,t') = T e (t') is allowed to depend on time t' . 
The situation with constant T e is well known from mean field spin glasses with one step 
of replica symmetry breaking PB[ , but we have not found such a constant T e beyond 
mean-field ||20||[2ll|. Only at exponential time-scales the mean field spin glass behaves 



as a realistic system 
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2.6. Time-scale arguments 

Consider a simple system that has only one type of processes (a processes), which falls 
out of equilibrium at some low T. When it ages a time t at T = it will have achieved 
a state with effective temperature T e , that can be estimated by equating time with the 
equilibrium time-scale, t = r eq {T e ). We have checked in solvable models that, to leading 
order in hit, it holds that T e = T e . (The first non-leading order turns out to be non- 
universal). This equality also is found in cooling trajectories, when the system is well 
inside the glassy regime. It says that the system basically has forgotten its history, and 
ages on its own, without caring about the actual temperature. Another way of saying 
is that dynamics in each new time-decade is basically independent of previous decade. 

This time-scale argument, however, is not very strong. Though it works in simple 
model glasses, it does not work, for instance, in realistic spin glasses. 



3. Solvable models 



In the above a variety of effective temperatures have been defined, and other definitions 
appeared also p71. The most prominent one seems to us the "thermodynamical" T e , 



that appears in the second law &Q < TdS cp + T e dX. If these effective temperatures are 
(basically) the same, then the above description leads to a coherent two-temperature 
picture. Let us now discuss models where this is or could be the case. (In principle also 
an effective field can occur [ lQj ; we shall not consider that complication.) 

The thermodynamic part of the above picture has been constructed from aging 
properties of the p-spin model ||24j| . If the limit iV — > oo is taken first, the aging 
dynamics starting from a random initial configuration follows from a set of equations 
very similar to the mode-coupling equations of glasses |24j. The long-time dynamical 
properties also follow from a replica calculation using the marginality condition [28[ . 
Our first step has been to find the physical meaning of the marginal replica free energy. 
Hereto we analyzed the TAP partition sum at zero field and non-zero field ]/]]. In 
both cases we convinced ourselves that the replica free energy (logarithm of the ordinary 
partition sum) indeed coincides with the logarithm of the TAP-partition sum. In other 
words, nothing went wrong by using replica's, and a physical interpretation should be 



anticipated. After the appearance of ref. |27] we were convinced that T/x should be 
interpreted as thermodynamic effective temperature |7|], and it was already known that 
this combination appeared in the fluctuation-dissipation relation |^SJ . Our studies with 
Sherrington and Hertz P0| on the p-spin model with a ferromagnetic coupling led to the 
insight that the fluctuation formula (|l6|) is valid there without the loss term. 

Similar results follow for a directed polymer model with glassy behavior f3~ll[p|1. In 
this model a directed polymer moves on a flat substrate with randomly located ridges. 
It prefers to lie in regions where ridges are widely separated. These Griffiths or Lifshitz 
singularities are the "TAP-states" of the problem. They occur often enough when the 
transversal width W scales exponentially in the longitudinal width W ~ exp(L 1 / 3 ). 
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The statics of this model has a Kauzmann transition. The long time dynamics involves 
hopping between these states. For an ensemble of non-interacting polymers starting 
from a random initial configuration, it has been argued that the motion involves a flow 
to states with decreasing complexity, as happens in the p-spin model at exponential 
timescales [[7]. Nevertheless, the dynamics of the model deserves further attention. 

Very informative are models with exactly solvable parallel Monte Carlo dynamics, 
such as independent harmonic oscillators |32| j^Of [^l|] or independent spherical spins in 
a random field p0|||2l|. In both cases the equilibrium timescale follows an Arrhenius 



law, which leads to glassy behavior at low temperatures when cooling or when aging 
from a random initial configuration. Since these models are not mean field-like, they 
have a proper dynamics, which satisfies the fluctuation formula ([16]). The effective 
temperatures decay as 1/lnt, with common prefactor, but different l/ln 2 t corrections. 

A model with a set of fast modes and a set of slow modes, that have a Vogel-Fulcher- 
Tammann-Hesse law for the divergence of the equilibrium time-scale, has also been 
formulated. It allows to study glassy dynamics below the Kauzmann temperature |10[ 



The two-temperature approach put forward here also explains that thermodynamics 
applies to black holes [[33| and star clusters |34| . 
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